Mode-Coupling Theory as a Mean-Field Description of the Glass Transition 
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Mode-coupling theory (MCT) is conjectured to be a mean-field description of dynamics of the 
structural glass transition and the replica theory to be its thermodynamic counterpart. However, 
the relationship between the two theories remains controversial and quantitative comparison is 
lacking. In this Letter, we investigate MCT for monatomic hard sphere fluids at arbitrary dimensions 
above three and compare the results with replica theory. We flnd grave discrepancies between the 
predictions of two theories. While MCT describes the nonergodic parameter quantitatively better 
than the replica theory in three dimension, it predicts a completely different dimension dependence 
of the dynamical transition point. We flnd it to be due to the pathological behavior of the nonergodic 
parameters derived from MCT, which exhibit negative tails in real space at high dimensions. 

PACS numbers: 64.70.qj, 61.43.Fs,64.70.pm, 66.30.hh 



The nature of the glass transition remains elusive de- 
spite of decades of discussion. Many theories and sce- 
narios have been proposed to explain the drastic slow 
dovi^n of dynamics of supercooled liquids but we still 
lack conclusive microscopic understanding of the phe- 
nomenon. Amongst various theories, mode-coupling the- 
ory (MCT) [ij and replica theory [2 Q are arguably the 
only first principles theories. 

On the one hand, MCT describes the slow dynamics 
of the mildly supercooled liquids using the static struc- 
ture factor as a sole input. It quantitatively captures the 
onset of the two-step relaxation of correlation functions, 
the scaling properties at the intermediate time scale (the 
P relaxation), and the algebraic increase of the struc- 
tural relaxation time. MCT, however, predicts a spuri- 
ous freezing transition at a lower density ipmct (or higher 
temperature Tmct) than the experimentally determined 
glass transition point (pg (or Tg). On the other hand, the 
replica theory is a static mean-field description of the 
glass transition '5]. It predicts that the fluid undergoes 
a thermodynamic or "ideal" glass transition at a higher 
density ipK than ipg (or lower temperature Tk than Tg), 
characterized by the one-step replica symmetry breaking. 
The replica theory also predicts that the dynamical tran- 
sition takes place at (pd < (pK (or Td > Tk) where the 
phase space or energy landscape starts splitting into nu- 
merous metastable states, or basins. MCT is conjectured 
to be the dynamical counterpart of the replica theory and 
V'mct to be identical to tpd, because mathematical struc- 
ture of MCT is equivalent to the dynamical equation of 
the so-called p-spin spherical model with p — 3, a. mean- 
field model for which the relation between the dynamical 
and ideal glass transition is rigorously established 0, Q . 
According to this mean-field scenario, the absence of the 
dynamic transition at (/Smct in real systems is interpreted 
as the round-off of the dynamic freezing by activated pro- 
cesses between basins in finite dimensions [sj . 

Despite the apparent and simple parallelism with spin 
glasses, the relationship of MCT with replica theory and 



physical insights from the mean-field treatments have 
never been fully understood. MCT was originally de- 
rived as a generalization of kinetic theories, using the pro- 
jection operator formalism with numerous uncontrolled 
approximations [IJ, whereas the replica theory is based 
on purely thermodynamic argument developed for disor- 
dered systems. In this Letter, we compare the MCT and 
replica theory results quantitatively, in order to clarify 
the relationship between the two theories developed in 
totally different arenas of physics communities. We es- 
pecially focus on the dimension dependence of the glass 
transition point and the nonergodic parameter fooiq), 
the plateau height of the density correlation function. 
To simplify the argument, we focus on the monatomic 
hard sphere system in d-dimension, for which the sole 
system parameter is the number density p — N/V or the 
volume fraction ip = Vdp, where Vd is the volume of a 
single hard sphere. We show that MCT is more quan- 
titative than the replica theory at c? = 3, which can be 
largely attributed to the lack of accurate approximation 
schemes in the replica theory to evaluate the static cor- 
relation functions of the replicated liquids. In higher di- 
mensions where the static correlation functions become 
trivial, discrepancies between the two theories become 
catastrophic. MCT's dynamical transition point {ip-mct) 
scales with dimension d differently from the replica coun- 
terpart ((/7d). This discrepancy comes from the spuri- 
ous negative tails of the van Hove correlation function, a 
generically positive quantity, that MCT predicts in high 
dimensions. This pathological negative tail is the origin 
of the non-Gaussian shape of /oo(<?) and thus the differ- 
ent d-dependence of (pmct from ipd- These results shed 
serious doubts over the validity of MCT in higher dimen- 
sions and call for reconsideration of MCT as a dynamic 
theory of the mean-field scenario of the glass transition. 

MCT is expressed as a set of nonlinear integro- 
differential equations for correlation functions such as 
the intermediate scattering function F{q,t) = 
{Sp{q,t)Sp{—q,0)), where Sp{q,t) is the density fluctu- 
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ation in reciprocal space at time t. The MCT equation 
for F{q,t) in d-dimension is given by 0,0] 

n-^F{q, t) + F{q, t)+ f ds M{q, t ~ s)F{q, s) - 0, (1) 



where Vlq = ksTq'^ /mS{q) is the phonon frequency 
and S{q) = F{q, i = 0) is the static structure factor. The 
memory function M{q,t) is given by 



Miq,t) 



dk 



\q+k\ 



dpV{q,k,p)F{k,t)F{p,t). (2) 



q — k\ 



In this expression, fc,p) ~ p S{q) s^-i J'^^^ kp 

{/+c(fc)+/_c(p)}^ /q'^+^ (47r)'' is the vertex function, 
where c{q) — p^^{l — 1/S{q)} is the direct correlation 
function, is the surface of d-dimensional unit sphere, 

is the Jacobian term. MCT predicts a nonergodic tran- 
sition at ipmct, beyond which the nonergodic parameter 
(NEP) /oo(g) = F{q,oo)/S{q) becomes non-zero. NEP 
can be calculated from the long time limit of Eq.([T|), 



/oo(g) 
l-/oo(<z) 



= Af(g,oo) 



(3) 



The self-part of the intermediate scattering function 
Fs{q,t) = {Sps{q,t)6ps{—q,0)) can be also described by 
an equation similar to Eq. ([T]) . The set of MCT equations 
can be solved numerically using S{q) as a sole input. 




FIG. 1: /oo(g) evaluated from MCT (solid line) and the 
replica theory (dashed line) at <^mct and ipd, respectively, q 
is scaled with the diameter of a sphere a. S{q) from HNC 
closure is used for MCT calculation. 

On the other hand, the replica liquid theory for the 
structural glass transition has been developed by Mezard 
and Parisi and recently applied to hard-sphere fluids 
by Parisi and Zamponi [3|. In this theory, replicated 
systems with a weak attractive interaction of order e be- 
tween them are considered. The free energy of the whole 
system is calculated as a function of intra- and inter- 
replica correlations. The relevant physical observables 



are evaluated by taking the limit e — at the end of 
the calculation. The theory predicts a thermodynamical 
transition characterized by the replica symmetry break- 
ing at which the configurational entropy vanishes. The 
transition point is often identified with the Kauzmann 
point ifK- The theory also predicts a dynamic transition 
point ifd < at which the free energy splits into numer- 
ous metastable basins. The system can not explore the 
whole phase space above (pd because of the infinite free- 
energy barrier which separates the basins. In the replica 
interpretation, the inter-replica pair density correlation 
function, g{r), plays the role of order parameter and is 
identified with NEP ph{q)/S{q) = foo{q), where h{q) is 
the wavevector representation of g(r) — 1. The NEP is de- 
termined from variational condition of the replicated free 
energy. In analogy with the p-spin spherical model [HI , it 
is believed that tpd should be identical to ipmct and that 
the NEP calculated from MCT should match that de- 
rived from the replica theory at the dynamic transition 
point. 

In order to check the validity of this conjecture, we nu- 
merically compare the results of both theories for (pmct, 
ipd, and foo{q) for d > 3. First, we look at a d = 3 hard 
sphere system. MCT for three-dimensional hard spheres 
was studied by Gotze et al. [1, On the other hand, 
the quantitative accuracy of the replica theory sensitively 
depends on the approximation scheme employed to calcu- 
late the free energy of replicated liquids. The small cage 
expansion technique is known to be a good approxima- 
tion near the Kauzmann point, but it does not describe 
the dynamic transition in low dimensions [H, [3| • There- 
fore, we use another scheme, the replicated hypernetted 
chain (RHNC) approximation the only method at 
present which captures the dynamic transition in finite 
dimensions. RHNC consists of a set of closure equations 
for both inter- and intra-pair density correlation func- 
tions, g{r) and g{r), given by 




ph{q) ' 
p[h{q) 



h{q)] 



p[h{q)-h{q)] J 
(4) 

Here, ph{q) — S{q) — 1 and v{r) is the interaction po- 
tential. The first equation is the HNC equation of a 
simple liquid [ll| and the second equation describes the 
inter-replica coupling. The dynamic transition point (pd 
is defined as the volume fraction beyond which h{q) 
becomes non-zero. We solve MCT equation, Eq.®, 
and the RHNC theory, Eq.Q, for the monatomic hard 
sphere system. For MCT calculation, we employ the 
HNC equation to evaluate S{q) in order to make the 
comparison consistent. The dynamic transition points 
thus obtained are Lp^act = 0.523 (which is slightly larger 
than 0.515 obtained from Percus-Yevick closure 9]) and 
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TABLE I: Values of (fmct and ipK from d = 4 to d = 8. ipK 
is from Parisi et al. 3]. 



d 


4 


5 


6 


7 


8 


fmct 


0.3652 
0.4319 


0.2542 
0.2894 


0.1736 
0.1883 


0.1159 
0.1194 


0.0751 
0.0739 



(pd ~ 0.612 [12|. fooil) calculated from MCT at ipmct and 
the replica theory at ifd are shown in Figure [TJ Quan- 
titative difference between the shape of the NEP from 
the two theories is obvious. Since it is well established 
that MCT's fooil) agrees very well with simulation 
and experimental results [l^, this discrepancy could be 
mainly due to poor performance of the replica theory. 
However, it is not clear whether this is attributed to 
the inherent inconsistency of MCT with replica theory 
or solely to a lack of accuracy of the RHNC approxima- 
tion. 

In order to give the two theories more stringent test, 
we discuss the dimension dependence of quantities near 
the dynamic transition point. We start with c? = 4 — 8 
and solve Eq.Q to evaluate (fimct and fooiq)- An algo- 
rism by Bans and Colot is used to evaluate S{q). 
<Pmct thus obtained is listed on Table U] along with (pK 
reported in Ref.[^]. Note that ipmct is smaller than tpx 
in lower dimensions, but the gap narrows with increasing 
dimension and at d = 8 ipmct exceeds ipK ■ 

Next, we study the d-dependence in even higher di- 
mensions, where the static properties of the liquid and 
the replicated liquid become insensitive to the approxi- 
mation schemes. Therefore, it is possible to check the re- 
lationship between the two theories without obscuration 
from approximations for the static inputs. In the high d- 
limit where the diagrammatic expansions of the free en- 
ergy is given by a simple function of the Mayer function 
Q-Pv(r) exact analytical expressions for static correla- 
tion functions are available. For the hard sphere system, 
the direct correlation function c{q) is simply given by 
c{q)=— {2'!T/qa)'''^^Jd/2iqcr), where Ji{x) is the Ith Bessel 
function of the first kind and a is the diameter of a sphere. 
Recently, accurate replica theory calculations of the free 
energy in high dimensions was carried for the monatomic 
hard sphere system, using this c{q) as an input and the 
cage expansion method [3| . The dynamic transition point 
was shown to scale with d as 



ipd = 4.8 X 2~''^d 



(5) 



and the Kauzmann point as (pK ~ 2^'^d\nd in the high 
dimension limit d — >■ oo Q . We solve the MCT equation 
with the same c(g) , keeping the convergence of discretiza- 
tion error and the numerical accuracy of the Bessel func- 
tion under control. In Figure [2l <^mct is shown as a func- 
tion of d. We find that (fimct scales as (/Jmct = 0.22x2^''(i^, 
in stark contrast with the replica prediction for ipd, 
Eq.©. We also calculate the NEPs from MCT and find 
that their g-dependence are non-Gaussian shaped in high 




FIG. 2: <^inct as a function of d. Filled circles are the nu- 
merical solution of Eq.(|3]) and empty circles from MCT with 
Gaussian approximation, Eq.(|6]). Solid line is the prediction 
from the replica theory. 

d. In low dimensions below d — 8, decay of ,foo{q) at 
qa > 10 and Js,oo{q) are well fitted by a Gaussian form, 
but in higher dimensions they decay faster than Gaus- 
sian at large g's (not shown). Different d-dependence of 
<y5mct from Lpd originates from this non-Gaussianity. This 
can be shown by solving Eq.Q assuming that foo{q) and 
.fs,oo{q) both have a Gaussian shape, i.e., foo{q), fs.,oo{q) 
« e~^« Z^''. Here, we assume that f^iq) « fs,oc{q), 
the so-called Vineyard approximation [11]. Substituting 
this Gaussian form in Eq.(|31), we obtain a self-consistent 
equation for R, 



1 
R 



psd 



2d^{2TTY Jq 



(6) 



This expression is strikingly analogous to the equation for 
the density field for the amorphous solid obtained using 
Gaussian approximation in the framework of the density 
functional theory |l5|. Solving Eq.([6]), we find that the 
equation has a finite solution for R above a volume frac- 
tion that behaves in high d as 



iflZl - 8.27 X 2-''d. 



(7) 



Though the prefactor differs, it retrieves the same d- 
dependence as the replica theory, which ascertains the 
origin of discrepancies between MCT and the replica the- 
ory. 

This non-Gaussian shape of the NEP that MCT pre- 
dicts in high dimensions is the compelling evidence that 
MCT breaks down. Analyzing the van Hove correla- 
tion function Gs{r,t), which is the real space represen- 
tation of Fs{q,t), makes this breakdown clear. Gs{r,t) 
is a distribution function of the distance for one par- 
ticle to explore during the time interval t, Gs(r,t) = 
{S{r - \R,{t) ~ R,{0)\)), where R,{t) is the position of 
the ith particle at time t. By definition, the van Hove 
function is a non- negative quantity. In Figure |3j we plot 
Gs.ooir) = Gs{r,t — oo) derived from MCT for several 
dimensions. Following the standard convention a mul- 
tiplicating factor Sdr'^~^ is used. In high dimensions, it 
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FIG. 3: The van Hove function Gs,oo{r) evaluated using MCT 
in several dimensions. From left to right, d — 4, 6, 8, 10, 15. 



exhibits a negative dip whose depth becomes larger as 
dimension increases. We checked that fs, ooil) retrieves 
the Gaussian shapes if the negative dips are absent. 

These results suggest that MCT in its present form 
is not consistent with replica theory and that moreover 
MCT suffers from serious deficiencies. The validity of 
the dimension dependence (fmct ^ that MCT pre- 

dicts is suspicious because it originates from the patho- 
logical behavior of the NEP. It is noteworthy that the 
non-Gaussian shape of fooil) and the negative tails of 
Gs.ooif) already appear ind — 8, below the upper-critical 
dimension dc = 8 of the glass transition 1^ 17 1 . Assum- 
ing a Gaussian shape for (q) in MCT recovers the lin- 
ear dimension dependence of 2 i^mct 

(see Eq.©) but the 
prefactor still does not match with the replica theory pre- 
diction. This fact implies that a quick remedy is unlikely 
to fix the problem. In hindsight, a convincing reason to 
conjecture that two theories are related is lacking, ex- 
cept for their apparent mathematical similarity with the 
p-spin spherical mean-field model of spin glasses. Decep- 
tively similar structure between Gaussian-approximated 
MCT, Eq.®, and the equation for the density profile 
derived in the mean-field analysis of the density func- 
tional theory [l^ also hints that MCT is a mean-field-like 
theory, but a small yet non-negligible difference between 
these equations leaves us a nagging suspicion over MCT's 
validity. One of routes to resolve these problems is to re- 
formulate MCT in a field theoretic languages in which 
parallels and differences with the mean field theory of 
spin glasses and the dynamic liquid theory is highlighted. 
Efforts in this direction have suffered from series of diffi- 
culties associated with consistencies with the fluctuation- 
dissipation theorem [l^l > the double-counting problem of 
the potential interactions, and the reconciliation between 
the dynamic and static liquid theories Although it 
is not clear if the replica theory and the density func- 
tional theory correctly describe the dynamic transition, 
at least our results clearly indicate that reconsideration 



and revision of MCT from the ground up are in order. We 
conjecture, however, that prospective revisions leave gen- 
eral mathematical properties of MCT equation intact. It 
is argued that MCT should be seen as a Landau theory 
in a sense that critical behavior and scaling properties 
that MCT describes near the dynamical transition point 
is universal [20j. Indeed, a recent numerical study in- 
dicates that MCT works better for critical behaviors in 
high dimensions [1]. 
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discussions. This work is supported by Grant-in- Aid for 
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Note added.- As this manuscript was being finalized for 
submission, we became aware of a paper by Schmid and 
Schilling [2l|. They have solved MCT equation for d > 10 
and shown the same dimension dependence of (pmct as 
reported in Figl2]and non-Gaussian shape of fooil)- 
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